Recoupling the non-secular part of nuclear spin-spin interaction in solids 
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NMR spin-spin relaxation in solids in strong magnetic fields is normally described only with 
the help of the secular part of the full spin-spin interaction Hamiltonian. This approximation is 
associated with the averaging of the spin-spin interaction over the fast motion of spins under the 
combined action of the static and the radio-frequency (rf) fields. Here we report a set of conditions 
(recouphng resonances) when the averaging over the above fast motion preserves some of the non- 
secular terms entering the full interaction Hamiltonian. These conditions relate the value of the 
static magnetic field with the frequency and the amplitude of the rf field. When the above conditions 
are satisfied, the effective recoupled spin-spin interaction Hamiltonian has an unconventional form 
with tunable parameters. These tunable Hamiltonian offers interesting possibilities to manipulate 
nuclear spins in solids and can shed new light on the fundamental properties of the nuclear spin-spin 
relaxation phenomenon. 

PACS numbers: 76.60.-k,76.60.Es 



I. INTRODUCTION 

The experimental technique of nuclear magnetic res- 
onance (NMR) has long established itself as a versatile 
tool for microscopic investigations of matter. One major 
source of microscopic information in NMR is the spin- 
spin relaxation. In solids, the description of NMR spin- 
spin relaxation requires averaging the full Hamiltonian of 
nuclear spin-spin interaction over the fast spin rotations 
induced by strong external magnetic fields. The present 
paper predicts theoretically a set of resonant conditions 
for solid-state NMR experiments leading to an unconven- 
tional form of the averaged interaction Hamiltonians. 

We consider nuclear spin-spin interactions in the pres- 
ence of a static magnetic field and a continuously applied 
rotating radio-frequency (rf ) field in the regime, when the 
two fields have comparable values, while both of them are 
still much larger than the local fields with which nuclear 
spins affect each other. Such a setting is realizable ex- 
perimentally. However, it implies the use of smaller-than- 
typical static fields and, therefore, a certain sacrifice in 
terms of the signal intensity. In a typical NMR exper- 
iment, the rf- field may be large in comparison with the 
local fields, but it is still much smaller than the static 
field. This "typical" setting for strong continuously ap- 
plied rf fields was originally considered by Redfieldi and 
subsequently discussed in NMR textbooks^i^. In such a 
case, the averaged interaction Hamiltonian is obtained 
from the full Hamiltonian with the help of a two-step 
truncation procedure^. The first step is to eliminate all 
interaction terms that do not commute with the Zeeman 
Hamiltonian for the static mag netic fieldi. The second 
step is to make the transformation into the rotating ref- 
erence frame where the rf field appears static, and then 
further truncate the interaction Hamiltonian by eliminat- 
ing all the terms that do not commute with the Zeeman 
Hamiltonian for the effective magnetic field in the rotat- 
ing frame. 

When the static and the rf fields are comparable to 



each other, the validity of the above two-step truncation 
becomes questionable, because the spin motions caused 
by both fields fall on the same timescale. The full inter- 
action Hamiltonian then needs to be averaged over the 
above two motions simultaneously rather than in two in- 
dependent steps. In this paper, we show that, in the case 
of a purely rotating rf field (as opposed to the case of 
a linearly-polarized rf field), such a single step averag- 
ing can be done rigorously without requiring the rf field 
to be smaller than the static field. We find that such a 
single-step procedure generically reproduces the results 
of Redfield's two-step truncation, but the reason for this 
agreement lies not in the perturbation theory but in the 
ergodic character of the single-spin dynamics driven by 
the static and the rf fields. We also find that in cer- 
tain special cases, which we call "recoupling resonances", 
the averaged spin-spin interaction Hamiltonian includes 
other terms in addition to the outcome of Redfield's trun- 
cation procedure. In the case of single spin species, the 
recoupling resonances appear when the rf field frequency 
is matched with certain multiples or simple fractions of 
the Larmor frequency corresponding to the effective mag- 
netic field in the rotating reference frame. In the case of 
two different spin species, the sums and the differences 
of the two effective Larmor frequencies also appear in 
some of the recoupling resonance conditions. The time- 
averaged interaction Hamiltonians associated with the 
above recoupling resonances are non-secular, in the sense 
that they do not conserve the total Zeeman energy asso- 
ciated with the external fields in either the laboratory or 
the rotating reference frame. These Hamiltonians may 
have a range of applications, such as, e.g., nuclear cross- 
polarization, or molecular structure determination. 

It is worth noting that the term "recoupling" is exten- 
sively used^"— in the context of the magic angle spinning 
(MAS) technique of NMR, where it implies the recovery 
of a part of the standard truncated Hamiltonian, after 
that Hamiltonian was suppressed by the spinning of the 
sample. The MAS recoupling schemes match the fre- 
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quency of the spin nutation^ or the frequency of the rf 
pulses^"— with the multiples of the sample spinning fre- 
quency. In the present paper, the term "recoupling" is 
used at a more basic level: it implies the recovery of a 
part of the full spin-spin interaction Hamiltonian lost in 
the standard truncation procedure. The recoupling reso- 
nances described below have nothing to do with sample 
spinning. However, they are best observable and, proba- 
bly, most useful in combination with MAS. 

In a broader context, our results could have been ob- 
tained in the framework of the averaged Hamiltonian 
theory^, but to the best of our knowledge it has never 
been attempted, perhaps because of the non-perturbative 
character of the problem. The role of the non-secular 
terms was, however, discussed in the context of the phe- 
nomenon of nuclear spin superradianceiSiii. 

The rest of this paper is organized as follows. We 
give the general formulation of the problem in Section [Hi 
introduce our single-step truncation procedure in Sec- 
tion mil present the results of this truncation with the 
identification of the recoupling resonances in Section IIVI 
and then discuss a possible experimental test and possible 
applications of the recoupling resonances in Section Ivl 



II. GENERAL FORMULATION 



of spins. The characteristic time of a single spin mo- 
tion induced by the Hamiltonian "Hd can be estimated 

/ \-l/2 

a.sT2 = hi J2j Jfj ) • The variable T2 is also to be re- 
ferred to as the transverse relaxation time. The timescale 
of nuclear spin relaxation due to the coupling to electrons 
or phonons is to be denoted as Ti (longitudinal relax- 
ation time). We further introduce the Larmor frequen- 
cies fioi = 7i-ffo a-iid the nutation frequencies uiu = 7ii?i. 
In the following, we drop lattice index i in the gyromag- 
netic ratio and use the notations 7 = 7,, i7o = 7-^0 1 
LJi = "/Hi, etc., when not dealing with the interactions 
of non-equivalent (unlike) spins explicitly. 

This paper is limited to the regime fioi'^i 3> I/T2 3> 
l/Ti. The first of these inequalities implies the separa- 
tion of the timescales between the fast motion due to the 
external fields and the slow motion due to the spin-spin 
interaction. The second inequality allows us to consider 
the nuclear spins as isolated from the environment. As 
mentioned in the introduction, the theoretical analysis 
in the literatureir— was limited so far by the additional 
inequality Hi <^ Hq (i.e. lui <C Oq). Our main focus 
is on the regime Hi ^ Hq. However, also, for a range 
of situations satisfying the inequality Hi <C Hq, the re- 
coupling resonances obtained below will imply important 
observable effects. 



We consider a lattice of nuclear spins in a static mag- 
netic field pointing along the z-axis, Hq = (0,0, Hq), 
irradiated by a rf field rotating in the x-y plane with 
frequency w, Hi{t) = {Hicos Lut, Hismu!t,0). We limit 
ourselves to the case of the magnetic dipole spin-spin in- 
teraction, which is, typically, dominant in solids. The full 
Hamiltonian^ii^ is 

n = n, + n^i + nD (i) 

where 

N 

n,^-HQY^"ia,„ (2) 

i 

N 

n,i = -Hi{t)-Y,"iJ,, (3) 

i 

and 

no = j:.h,[i.-i,-^-^^^^^^^^^]. (4) 

Here J;,,- — — are the coupling constants, li = 

{lixi liyi liz) are the operators of the nuclear spin of the 
ith lattice site, 7^ is the gyromagnetic ratio of that spin, 
Tij = {rij^x,Tij^y,rij^z) are the displacement vectors be- 
tween two lattice sites i and j, and N is the number 



III. TRUNCATION PROCEDURE 

We calculate the effective spin-spin interaction Hamil- 
tonian by averaging the full interaction Hamiltonian % d 
over the fast spin motion induced by the terms T-Lz +'Hj.f. 
In the Heisenberg representation, this fast motion is de- 
scribed by the following equation for each of the spins : 

^=/,X7H(0, (5) 

where H{t) = Hq + Hi{t). The above equation is lin- 
ear in terms of spin operators {lix, hy, liz), and, there- 
fore, its predictions are essentially the same as those for 
the classical spins, in which case {lix, hy, liz) would be 
simply three numbers. Despite the linearity of Eq.([S]), 
it is not solvable analytically either classically or quan- 
tum mechanically for an arbitrary time-dependent H{t). 
However, for the specific time dependence considered in 
this work, the solution is facilitated by the transforma- 
tion from the laboratory reference frame into the refer- 
ence frame rotating around the z-axis with frequency lo. 
The equations of motion in the rotating frame have the 
same form as Eq. ([5]) , but with time-independent effective 
magnetic field He = {Hi,0,Hq +uj/"f). (The x-axis in 
the rotating frame coincides with the direction of the rf 
field.) The corresponding Larmor frequency is 

rie = 7i/e = sign(7) ^ujI + (flQ+Uj)^. (6) 
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Thus, the motion of a spin in the laboratory frame can 
be characterized as the Larmor precession with frequency 
He around the direction of He, which itself rotates with 
frequency w around the z-axis. 

Our formal procedure for averaging Ji d consists of the 
transformation into the reference frame that follows the 
fast spin rotation induced by Eq.(I5]). In that reference 
frame, each spin would be motionless, if the full Hamil- 
tonian included only Ji^ + "Hj-f . Such a transformation 
is time-dependent. As a result, the Hamiltonian "Hd ex- 
pressed in the new spin coordinates acquires some quickly 
oscillating interaction coefficients, which we then average 
out. 

The above procedure may, at first sight, appear equiv- 
alent to the one introduced by Redfield^"'^. The difl^er- 
ence, however, is that Redfield did not define the entire 
transformation and then average the resulting Hamil- 
tonian over time but rather used the Zeeman-energy- 
conservation criterion of Van Vlecki to average sepa- 
rately the Hamiltonians obtained at the intermediate 
steps of this transformation. Such an independent av- 
eraging of the intermediate Hamiltonians is justifiable by 
the perturbation theory argumentsi only when f^o 3> Wi . 
It is shown below, that in the case of CIq ^ uji, the spin 
motions associated with the intermediate transformation 
steps may be correlated, and thus the averaging over 
them cannot be done independently. 

The overall transformation into the reference frame 
that follows the motion induced by Eq.(l5]) can be de- 
composed into three simple transformations visualized in 
FiglH The first transformation is from the laboratory 
frame to the reference frame rotating with the rf field: 



cos(cji) sin(cji) 0"* 
Uu, = I — sin(cji) cos(cji) 
1, 



(7) 



The second transformation is from the rotating to the 
"tilted-rotating" reference frame, where the z-axis coin- 
cides with the direction of He'. 



'cos a —sin a' 
C/a I 1 
, sin a cos a 



with 



a = arccot 



no 



(8) 



(9) 



admitting values in the range [0,7r]. 

Finally, the third transformation is into the "double- 
rotating" reference frame precessing around Hf, with fre- 
quency — fie'. 



^coaiflet) -siniflet) 0^ 
Uq^ = \ sm{flet) cos(f2ei) 
1, 



(10) 




Figure 1: (a) Illustration of the transformation from the 
laboratory reference frame (a;, y, z) to the rotating reference 
frame {x^ , z^) — Eq.([7]). (b) Illustration for the transfor- 
mation from the rotating reference frame (x"^ , z^) to the 
tilted-rotating and then double-rotating reference frames — 
Eqs.dHHO}. 



The effective magnetic field in this reference frame is 
equal to zero. The explicit implementation of trans- 
formation (|10p before performing the time-averaging of 
Hd formally discriminates our treatment from that of 
Redfieldl-^. 

We denote the axes in the double-rotated reference 
frame as {a;', y', z'}, and the corresponding spin projec- 
tions as II = {lix, liy, I'iz)- The spin coordinates in the 
laboratory reference frame can now be expressed as 



(11) 



The substitution of Eq. (fTT|) into Eq. @ allows us to reex- 
press the interaction Hamiltonian % d hi terms of I[ and 
then average the interaction coefficients in the resulting 
expression. 

The full time-dependent form oi %£> in the double- 
rotated reference frame is rather long. It can be found 
in Ref together with the implementation of the time- 
averaging procedure. Below we exemplify this procedure 
by one typical term and then, in the next section, present 
only the resulting time-averaged expressions for T-Ld- 

The term we have chosen as an example, originates 
from the calculation for equivalent (like) spin species dis- 
cussed in the next section. It has form 



j{t) c{rij) Ij^^L^, 



(12) 



where c{rij ) is some function of the relative displacement 
between the two nuclei, and 



aij{t) = cos^ {(jjt) cos^ (flet) ■ 



(13) 



The time-averaged value of the above coefficient is calcu- 
lated as follows 



= lim 



aij{t)dt 



i ifw^±r!e, 
I ].i0J = ±Vte. 



(14) 




In the above case, (oy) = |; is what one obtains by 
the standard two-step truncation procedure outhned in 
the introduction, while, in the case of w = if^e, there 
is an extra contribution associated with the fact that 

(C0s2(wt)c0s2(f}et)> ^ (C0s2(wt))(c0s2(f}et)). 

The question then arises, what happens, if, in an ex- 
periment, the actual values of _ffo, Hi and w are such 
that one of the recoupling resonance conditions is satis- 
fied approximately but not exactly. In such a case, the 
Hamiltonian-averaging procedure needs to be modified. 
Below, we illustrate this modification for the resonant 
condition lo = ^\^e but the expression for the resulting 
correction to the averaged Hamiltonian is applicable to 
all other recoupling resonances. 

In the chosen example, the mismatch of the resonant 
condition can be parameterized with the help of vari- 
able Afig defined by equation uj — — -f- Ailg). Now 
we modify the last step of the transformation to the 
double- rotated reference frame [Eq. (|10p ] by changing the 
value of the rotation frequency from f^g given by Eq. ^ 
to + Afig. In the modified double-rotated reference 
frame, the averaged interaction Hamiltonian obtained in 
Section IIV Al for the recoupling resonance a; = 
will be exactly valid, but, in addition, since f2g + Af2g 
is different from Jig, the effective magnetic field will be 
different from zero, which then gives rise to the following 
additional term: 



A?^ = Af7g^4. 



(15) 



When Ailg ^ 1/72, such a term suppresses the ef- 
fect of the resonant non-secular terms given in the next 
section. This is, of course, expected for a significant de- 
parture from the recoupling resonance conditions. When 
Arig ^ ^IT2, the above term still suppresses the effect 
of the non-secular terms but only partially, and this par- 
tial suppression is of the same order of magnitude as the 
suppression caused by the secular interaction terms (to 
be explained in Section FlV A[) . In such a case, the term 
(jlSp should simply be included in the averaged Hamil- 
tonian together with the secular terms and the resonant 
non-secular terms. 

The above considerations imply, that, in a generic case, 
when the strength of the secular terms can be estimated 
as h/T2, the recoupling resonance conditions cannot be 
resolved experimentally with accuracy much better than 
±1/T2. However, if the secular terms are suppressed by 
other means, then the experimental resolution will be 
limited by the strength of the recoupled terms them- 
selves. Thus, if the strength of the recoupled terms is 
much smaller than S/T2, then the recoupling resonances 
can be resolvable with proportionally better accuracy. 

Let us now think about the averaging procedure in 
terms of classical spin trajectories. The typical situa- 
tion corresponds to an irrational value of the ratio w/Jlg. 
In this case, the trajectory of a classical spin that fol- 
lows Eq.([5]) is open and ergodic on the spin sphere, as 
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(d) 



Figure 2: Examples of trajectories of a single classical spin 
under the combined action of the static and tiie rf magnetic 
fields: (a) general case, open ergodic trajectory; (b) uj = — f2e, 
closed trajectory; (c) u) = —2Q,e, closed trajectory; (d) uj = 
if^e, closed trajectory. 



illustrated in Fig. Ufa). Averaging over such trajecto- 
ries leads to the standard two-step truncation result. 
If w/rig is a rational number, then the individual spin 
trajectories become closed, and one can suspect that 
the averaging result is different from the irrational case. 
In reality, however, almost all rational values of w/fig 
lead to the same averaging results as the irrational val- 
ues. In the case of like spins, the only exceptions are 
w/rjg = {±1/2,-1,-2}— see Figs.d^b-d). The fact that 
no other rational numbers appear in this set is related to 
the fact that the interaction Hamiltonian T-Ld is of the 
second order in terms of spin variables. If higher order 
spin couplings were involved, e.g. of the form hxljxlky, 
then other ratios with absolute values such as 1/3 or 3 
may also appear in the above set. 



IV. RESULTS 



In this section, we present the time-averaged Hamil- 
tonians separately for the couplings of "like" spins (sub- 
section IIV A[) and "unlike" spins (subsection IIV Bp . Like 
spins have the same gyromagnetic ratios and hence the 
same Larnior frequencies. Unlike spins have different gy- 
romagnetic ratios and different Larmor frequencies. The 
Hamiltonians referring to the above two cases are marked 
by superscripts ' and respectively. The details of the 
calculations can be found in Ref J^. 
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A. Like spins 

Away from the recoupling resonances, our averaging 
procedure reproduces the standard truncated Hamiho- 
nian for hke spins^: 



DO 



1 



2 

COS a 



I' I' 



^{lUj.+nyljy) , (16) 



where a is given by Eq. ([9]) . 

As aheady mentioned in Section IIIIl the averaged 
Hamiltonian contains additional couphng terms in the 
case of recouphng resonances associated with ratios 
u/^e = {±1/2,-1,-2}. According to Eq.®, these ra- 
tios translate into the following expressions for u in terms 
of Qq and uJi: 



2^0 



U} = ± — J7e 

2 



(17) 
(18) 
(19) 



The upper (lower) sign before the square root in Eq. (|19p 
corresponds to the upper (lower) sign of the resonant 
condition u; = i^ile. In the case of Eq. (fT5|) both signs 
before the square root are realizable for the same con- 



dition 



-2i7e as long as 



> 



V3\ 



UJl\ 



When the 



value of uj following from Eg . P?|) . (fT5|) or (|19p is negative, 
this indicates that frequency vector u: is anti-parallel to 
Hq. The resonance conditions to = fig and w = 2Qe are 
absent in the above list, because, after the substitution 
of given by Eq. ^ they give no real- valued solutions 
for Lj. 



Now we list the corresponding averaged Hamiltonians: 



-fie 



- ( COS Q — 1 ) 



-il' I' 



I' I' ) 



+ 



2 '''ij,x^ij.y 



(I' I' + I' I' )] 



+ -(cos 2a ■ 



cos a) 



'-d' I' 

y-'-ix-'-j 



+ 1' I' ) 



:fj' /'. +/' J' -ill- 



(20) 



to — — 2f2e 



'T-D2 — rtoo 



5^ Jij — [sin 2Qf — 2 sin a] 



I' _ J' J' ) 

y-'-ix-'-jx ^y jy J 



LU = ± — f2e 



ri-D3 — rLoo 



[sin 2a ± 2 sin a] 



Itr' /'. +1' I' )^l ^ ^ij,^^^j,y 

y-'-ix-'-jz ^ -'-tz-'-jxJ 2 



(4/;. + 44)]- (22) 



The upper (lower) sign of the resonance condition of 
Eq. (|22p corresponds to the upper (lower) sign in the ex- 
pression for ^^3. 

We note that the standard truncated Hamiltonian 
conserves the z'-projection of the total spin polarization 
J2i I'iz tl^6 double-rotating reference frame. This to- 
tal spin component still relaxes on the time scale of the 
order of Ti because of the interaction with electronic 
spins or phonons. We call this "normal longitudinal re- 
laxation." In contrast, the recoupled terms appearing in 
Eqs. pOII22|) do not conserve Z-^, which means that all 
three projections of the total spin polarization decay to 



zero in both the double-rotating and the laboratory ref- 
erence frames on a time scale, which is, in general, much 
faster than Ti. We call this "anomalous longitudinal re- 
laxation." 

The recoupled terms in Eqs. (|20j|2^ are comparable 
with H^g, when Hi ^ Hq. In such a case, the time scale 
of the anomalous longitudinal relaxation is expected to 
be of the order of T2, which is clearly much faster than 
Ti . Therefore, the recoupling resonances should be read- 
ily identifiable in an experiment of the type proposed in 
Section rVAl 

In the case of Hi <C iZo , the recoupled terms can still 



6 



induce the anomalously fast longitudinal relaxation, but 
this case requires a more detailed discussion. We first 
note that all recoupled terms in Eqs. (|20ll2^ are equal 
to zero, when a = 0, which, according to Eq.Q, cor- 
responds to Hi ~ 0. When Hi ^ Hq, the recoupled 
terms are much smaller than the terms in "H^q, in agree- 
ment with the standard truncation result of Redfieldi"— . 
The values of the a-dependent prefactors of the recou- 
pled terms in the Hamiltonians (|20j|22p are plotted as 
functions of Hi/Hq in FigO for Hi/Hq < 0.5. The 
strongest of these terms are of the order of Hi / Hq . They 
correspond to the resonances lo ~ and lu — 257o 
obtained, respectively, from Eg. ljlSp with plus sign and 
Eq. (|19p with minus sign in the limit Hi/Hq — > 0. (In om 
sign convention, the usual single-spin NMR resonance is 
located at u! « — ^io-) 

In general, the recoupled terms linear in Hi/Hq should 
induce a longitudinal relaxation on the time scale of the 
order of T2{Hq/Hi)'^, which can still be much shorter 
than Ti and thus represent a clear indication of a re- 
coupling resonance. The above estimate is of the second 
order in terms of Hq/Hi, because the first-order effect 
of the recoupled terms should be averaged out by the 
spin motions due to the much stronger secular terms. In 
view of this consideration, the observability of the recou- 
pling resonances can be significantly improved by sup- 
pressing the secular Hamiltonian H^q. If this is done, 
the timescale of the longitudinal relaxation induced by 
the recoupled terms can be estimated as T2{Hq/ Hi), i.e. 
it becomes much faster. 

The suppression of "H^q makes the observable effects 
of the recoupling resonances more dramatic not only in 
the limit of Hi ^ Hq, but also in the general case of 
Hi ^ Hq. It can be realized, for example, as follows. 

A given recoupling resonance imposes a constraint on 
the three parameters Hq, Hi and uj, which leaves a free- 
dom to impose one additional relationship. This free- 
dom can be used to impose the magic angle condition 
1 — 3 cos^ a — Q. However, this condition together with 
any of the like-spin recoupling resonances fixes the ratio 
Hi/Hq to be of the order of 1, which means that it can- 
not be used if the experiment is limited to the regime 

Hi < Hq. 

Another resource for suppressing is to exploit the 

spatially dependent factor 1 For example, if the 

recoupling resonances are to be investigated in CaF2, 
where ^^F nuclei form a simple cubic lattice, it is bet- 
ter to use the static field along the [111] direction of that 

lattice. In such a case case, 1 = for all six near- 

est neighbors of each spin, i.e. the secular interaction 
with these neighbors is suppressed, while the recoupled 
terms remain large. 

Of course, the most effective way to suppress the secu- 
lar interactions for all pairs of interacting spins is to use 
the MAS technique^. This technique completely aver- 

ages out all factors 1 :/^, when cos^ Q = 1/3, where 




3/4[sin(2»)-2sin(ff)| w+=-2ne ~ -2/3no 

\ 3/4|sin(2ff)-2sin{»)| 

^ \ , , ; \ — L 



0.001 0.005 0.010 0.050 0.100 0.500 

H^|Ha 

Figure 3: (Color online) Plots of the a-dependent prefactors 
of the recoupled terms in the like-spin Hamiltonians (|20I22|) as 
functions of the ratio Hi/Hq for Hi < 0.5Ho. The legend for 
each line indicates the corresponding resonant condition and 
the a-dependence of the prefactor being plotted. Each reso- 
nant condition is given both in its exact form (|17I19[) and in 
the approximate form for the case Hi <C Hq. The conditions 
LJ+ — — 2f2e and = —2Qs correspond to the recoupling 
resonances (|18l) with signs (+) or ( — ), respectively. 



Q is the "magic" angle between the spinning axis and 
the static magnetic field. Spinning at the magic an- 
gle, however, does not average out completely the factors 
^ij.y-'^fj,xi '^ij,xnj,y, Or rij^^Tij^^ that enter the recoupled 
Hamiltonians. 



B. Unlike spins 

In this subsection, we consider coupling between two 
unlike spin species. We will use lattice index i to label 
one of these species and index j to label the other. We 
reintroduce the labeling of gyromagnetic ratios 7^ and 7^ , 
and the corresponding Larmor frequencies VLqi — jiHQ 
and Oqj ~ IjHq. The two spin species are now to be de- 
scribed in two different double-rotated reference frames. 
The effective Larmor frequencies ^lei, ^ej and the corre- 
sponding tilting angles ai , aj are then defined by adding 
the indices i or j to variable f2e, f^o, uji and a in Eqs.® 
and ([H). 

As in the case of like spins, our averaging proce- 
dure for unlike spins recovers the standard truncated 
Hamiltonian^ away from the recoupling resonances: 

^L>o ~ Jij-jl^cosai costtj — cos(ai — aj)] 

X (1-^)4/;., (23) 

where the sum j) includes all pairs of unlike spins. 
It should be mentioned here, that, in the very special 
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case of 7i = — 7^ and w = 0, the secular Hamiltonian 
T-l^Q should also include the so-called "double-flip" terms 
and Z-./j-, where I'^^ = I^^ ± I[y, etc. These 
terms are preserved by the averaging, because their coun- 
terparts li^Ij^ and li-Ij- in the original full Hamil- 
tonian ([T]) now conserve the Zeeman energy. This case 
can be mapped onto a somewhat unusual problem of like 
spins by changing the sign of all y- and z- spin projec- 

~Ijy and 



tions for one of the two spin species, e.g. Ij 



Ijz — >■ —Ijz- However, when w 7^ 0, the two spin species 
in the transformed problem will experience rf-fields ro- 
tating in the opposite directions, which means that the 
transformed problem is not equivalent to the problem of 
like spins in a single rotating field. 

Another special case is 7^ — )• 0. It implies that the 
terms of the type Iizlj± in the full Hamiltonian ([1]) be- 
come secular, and, as a result, terms Iizlj± may appear 
in 'H^Q. Such a limit is not realizable physically in the 
case of magnetic dipole interaction, because 7^ = si- 
multaneously suppresses the interaction itself. 

The above special cases should be kept in mind in the 
analysis of the parameter dependence of the recoupling 
resonances. 

There exist two kinds of "unlike" recoupling resonances, 
namely, the resonances involving the motion of only one 
spin species, and the resonances involving both spin 
species. The resonances of the first kind are the following: 



UJ = ± — fip 

2 



(24) 
(25) 



and also those obtained from Ea. (|24l25p by replacing in- 
dex i with index j. The resonances of the second kind 
are 



1 



(-l)"-(l)e.±r!e,), 

-1) 



(-ly 



2 

'in. 



(26) 
(27) 



where variable n takes values or 1 and thus controls 
the sign in front of the entire expression. Resonant con- 
ditions (|26l27p lead to fourth-order polynomial equations 
with respect to w, which, in principle can be solved an- 
alytically. However, in view of the cumbersome char- 
acter of the resulting formulas, we chose to investigate 
the solutions numerically. Several plots illustrating the 
character of these solutions are presented in Fig. |4] and 
further in Fig. [5] of Appendix [Appendix A[ As evident 
from these plots, it cannot be guaranteed that each of 
the eight conditions envisioned by Eas. (|26l27p is realiz- 
able for a given set of values of 7^ , 7^ , Hq and Hi . For 
example, if 7^,7^ > 0, then no value of uj exists that 
would satisfy the condition uj = l/2{ilei + ilej)- 



The time-averaged Hamiltonians corresponding to the 
recoupling resonances (|26l27p are the following: 



I/" — 11^ 



,j- [cos{ai + aj) 



2 ^ix-'jz "T 



I' I' 1 



(28) 



UJ — ±-Qei 



riD: 



H 



'J.a: r' t' _ 



2 r, 



ij,x' ij,y ji ji 



(29) 



UJ = (— l)"(r2ei ± Qej) 



"Hm = y-DO + X] Jij^ [sin(ai -I- aj) ± (-l)"sinaj + (-l)"sinQj] 



(•/' J' =F/' /' \ + (±1' /' +7' /' \ 



(30) 



UJ = (-l)"-(nei ± Qej) 



rtn^ 



"Hdo + [cOSQi -I- (-1)"] [cOSOj ± (-1)"] 



'lU-xTliyl'jy)-{-l) 



,2r. 



n ^ ' ij,xi tj.y 



±J' J' + /' J' ) 



(31) 



The upper (lower) sign of the recoupling resonances of Eqs. p0|3ip corresponds to the upper (lower) signs in the 
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Figure 4: (Color online) Examples of recoupling resonances 
for unlike spins. Plots of the resonant values of uj and the 
prefactors in the recoupled Hamiltonians pOI31[l as functions 
of 7j/7i for two resonant conditions indicated above the re- 
spective plots. All examples presented are obtained numeri- 
cally for either Hi/Hq = 0.9 or Hi/Hq ~ 0.1. Each column 
of plots corresponds to the ratio Hi /Ho indicated above it. 
The plot lines consist of the actual computed points and thus 
in some cases have the appearance of dotted lines. The re- 
sults for Hi /Ho = 0.1 are qualitatively representative of the 
limit Hi <C Ho- The analogous plots for the six remaining 
resonant conditions associated with Eas. (|30l3ip are given in 
Appendix [Appendix A| 



for at least one of the two spin species. We observe, 
however, that the total z'-polarization of the both spin 
species, i.e. J2i ^iz + J2j ^jz^ is still conserved by the re- 
coupled Hamiltonians given by Ea. ([30| for u = zt{ilei — 
ilej) and by Ea. ([3T|) for uj = ±-^{flei — ^ej), because the 
terms in these Hamiltonians have the character of the 
so-called "flip-flops." 

In general, the recoupled Hamiltonians ()24j|27p are 
comparable to the standard truncated Hamiltonian , 
when Hi is comparable to Hq — see plots in Fig. 3] for 
Hi/Hq — 0.9. If Hi <C Hq, then the recoupled terms 
are normally small, because Eq.([9]) gives the values of 
ai and aj either close to or to tt — Figs. [5ljc,d,f) for 
Hi/Ho = 0.1. 

However, as the plots in Figs. |4] and [5Ua,b,e) for 
Hi/Hq = 0.1 indicate, and our approximate calcula- 
tions confirm, there are interesting and potentially use- 
ful exceptions, when the angle-dependent prefactors of 
the recoupled Hamiltonians remain comparable to 1 even 
though Hi <C Hq. The first exception corresponds to 
the condition Jj/ji ~ —1 for the recoupling resonances 
UJ = ±l/2(ilei — ^lej) and uj = l/2{ilei + ^ej)- In these 
cases, the recoupled terms in Hamiltonians pip remain 
large as long as |7j/7j + 1| < Hi/Hq. [In the case of 
resonance u ~ l/2{flei + ^ej), one should also be mind- 
ful of adding the secular double-flip contribution men- 
tioned after Eq. (|23|) .] This group of resonances can thus 
be exploited for pairs of nuclei satisfying the condition 
7i w — 7j, such as e.g. ^■^^Xe and ^^Na. The second ex- 
ception corresponds to the condition of Jj/ji ~ (or 
7j/7j w 0) for the resonances lu = ±{^ei + ^ej)- The 
corresponding recoupled Hamiltonians (|30p are large, as 
long as |7j/7i| < Hi/Hq. These recoupling resonances 
can be very useful for the purposes of cross-polarization, 
when the gyromagnetic ratio of one of the two nuclei is 
much smaller than the other. 

Finally, we mention that, as in the case of like spins, 
the suppression of the secular interaction terms between 
unlike spins by matching the magic angle condition or by 
the MAS technique is an additional resource for observing 
the recoupling resonances. The only difference in the 
present case is that the magic angle condition involves 
two angles: Scoscti cosaj — cos(Q;i — ctj) = 0. 



V. DISCUSSION 



A. Experimental realizability 



corresponding Hamiltonians. In the case of a small mis- 
match of the above resonant conditions, the additional 
term psp for one or both spin species should be included. 
In the case of resonances (|30l3ip . one, in fact, has the 
freedom of distributing the mismatch between the two 
spin species. 

The standard truncated Hamiltonian conserves 
both ^I'z '^j^jz^ while each of the recoupled 
Hamiltonians (P5)) - ([3T|) violates the above conservation 



The standard truncated Hamiltonians (|16l23p conserve 
the z'-projection of the total spin polarization for each 
of the available spin species in their respective double- 
rotated reference frames. The most obvious qualita- 
tive effect of the recoupled Hamiltonians is that they do 
not conserve this total z'-projection for at least one of 
the spin species. This difference can be tested experi- 
mentally. In the case of like spins, it will lead to the 
anomalously fast longitudinal relaxation discussed in Sec- 
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tion llV Al In the case of unlike spins, it may, in addition, 
lead to an anomalous cross-polarization effect. 

To be specific, a possible way to observe the anomalous 
longitudinal relaxation can consist of the following steps: 
(i) polarizing the system in the static magnetic field; (ii) 
turning on the rf field (preferably, adiabatically) to match 
one of the recoupling resonance conditions; (iii) waiting 
enough time for the additional interaction terms to relax 
J2i I'iz (the waiting time should still be much shorter than 
Ti); and, finally, (iv) suddenly turning off the rf field. No 
free induction decay signal should be observable in this 
case. In contrast, away from the recoupling resonances, 
a free induction decay signal of finite intensity should 
always be observable for waiting times of the order of Ti . 

The requirements for the above kind of experiments 
are the following: (1) The static and the rf fields should 
be strong enough to ensure VLq.uji ^ I/T2. (2) As 
discussed in the context of Eg. ljlSp . the recoupling res- 
onances should be matched to an accuracy of roughly 
(3) The static and the rf magnetic fields should be 
relatively homogeneous; otherwise the recoupling condi- 
tions cannot be satisfied over the entire sample. (4) The 
strength of the recoupled terms must be such that the 
timescale associated with the anomalous relaxation of the 
longitudinal magnetization is faster than Ti (see the end 
of Section flV A[) . (5) Finally, in the case of Hi ^ Hq, a 
real circularly polarized rf field should be used, as done, 
e.g., in Refji^. 

Now we discuss the observability of the recoupling res- 
onances in the NMR setting that uses a linearly polarized 
rf field in the regime of Hi <C Hq. The latter condition 
extends to the experiments in larger static fields, which 
means stronger NMR signals. 

The Hamiltonian averaging procedure used in this pa- 
per cannot be straightforwardly extended to linearly po- 
larized rf fields, because a closed analytical solution for 
the spin motion is not available in such a case. In par- 
ticular, the usual assumption that one can neglect one of 
the two counter-rotating components contributing to the 
linearly polarized rf field is not justified when Hi ^ Hq. 
However, when Hi <g; Hq, the case of linearly polarized 
rf field should be treatable by combining our truncation 
procedure with the perturbation expansioii^iii^. In this 
limit, the presence of two rotating components may lead 
to additional recoupling resonances, but the resonances 
obtained in the present paper for a single rotating com- 
ponent of the same frequency should also be present. We 
further expect that the same recoupling resonance in the 
rotating and the linearly polarized settings is generically 
characterized by the same recoupled Hamiltonian, when 
the leading order of the expansion of this Hamiltonian in 
terms of Hi/Hq is linear. 

At small values of Hi/Hq, the recoupling resonances 
satisfying the above linearity requirement are also the 
strongest and thus the most promising in terms of ex- 
perimental observation. As discussed in Section IIV Al for 
the case of like spins, the linearity with respect to small 
Hi /Hq is exhibited by some but not all recoupled Hamil- 



tonians. We also note in this context, that an experiment 
with small linearly polarized fields should avoid the inter- 
play with the secondary single-spin resonances described 
by Winter^ii^.) 

In addition to the possibility of using stronger static 
fields and linearly polarized rf fields, the regime of Hi <t; 
Hq has two other advantages. First, the rf field homo- 
geneity requirement is less stringent, when Hi <C Hq, 
because what matters is the homogeneity of the effective 
frequencies Jleij to which the leading rf field contribu- 
tion is only of the order of (Hi/Hq)'^. Second, in the 
case, when the secular part of the interaction is exter- 
nally suppressed, the smallness of Hi/Hq implies that 
the recoupling resonances can be resolved with higher 
accuracy (see the discussion in Section HlH) . 



B. Possible applications 

Fundamental studies of spin-spin relaxation: Perhaps, 
the most direct use of the tunable Hamiltonians (|20II221 
I28ll3ip is to conduct fundamental experimental studies of 
nuclear spin-spin relaxation in solids, which is still not 
completely understood, in a much broader range of pa- 
rameters than those accessible with the standard trun- 
cated Hamiltonians. 

Cross-polarization: In the case of unlike spins, the re- 
coupled Hamiltonians (|30ll3ip for the resonances uj = 
±(f2ei — ftej) a-nd = ^^(^ei ~ ^ej) Can be used to 
cross-polarize different spin species. As mentioned in Sec- 
tion IIV Bl these Hamiltonians contain "flip-flop" terms 

^'iy'^jy ^hich can transfer 
polarization from one spin species to the other. Such a 
transfer would require an rf fleld rotating with a single 
frequency, as opposed to the standard Hartmann-Halm 
cross polarization routinei^, which involves two rf-field 
frequencies. The two spin species can also be cross- 
polarized, albeit less efficiently, using the transient effect 
of the "double-flip" terms I'iJ'jy+Iiyl'j.^ and I[J[j^~I[yI[jy 
for resonances lo = i{flei + ^ej) a-nd u = + ^ej)- 

The suppression of the secular terms by either the magic 
angle condition or the MAS technique, as discussed in 
section IIVBI can enhance the efficiency of both the flip- 
fiop-based and the double-fiip-based cross-polarization 
processes. In the latter case this enhancement should 
be particularly significant. 

Structure determination: The use of solid-state NMR 
for determining the structures of complex molecules in- 
volves MAS decoupling between nuclear spins followed by 
selective recouplin g^^i^^ . The recoupling resonances ob- 
tained in this paper can, possibly, be used at the latter 
stage. As discussed in Section |VA| the regime Hi <C Hq 
may provide both the sufficient intensity and the suffi- 
cient frequency resolution for such experiments. If suc- 
cessful, such a method will not require the recoupling se- 
quence to be adjusted to the sample spinning frequency. 
It will also have an added benefit of the non-trivial spa- 
tial dependence in the recoupled Hamiltonians as com- 
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pared to the usual secular interactions. This dependence 
can be converted into additional structure information 
for oriented samples. 

In summary, the recoupling resonances for like and un- 
like spins amount to a new resource for manipulating 
nuclear spins, and, therefore, are likely to be useful for 
various NMR applications. However, the condition of 
applying a high amplitude rf field for a rather long time 
may be challenging to realize in practice. 



VI. CONCLUSION 

In this paper, we described theoretically the properties 
of recoupling resonances associated with the appearance 
of non-secular terms in the time-averaged Hamiltonians 
of nuclear spins-spin interaction in solids in the presence 
of strong rotating rf fields. Our analysis indicates that, 
under certain conditions, the effect of the recoupled non- 
secular interaction terms should be particularly well ob- 



servable. We have also discussed the extension of our 
findings to the case of linearly polarized rf fields. Finally, 
we discussed how the recoupling resonances considered 
in this paper can used for manipulating nuclear spins in 
various NMR applications. 
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Appendix A. PLOTS OF RECOUPLING 
RESONANCES FOR UNLIKE SPINS 

This appendix contains the collection of plots (Fig. [5]) 
for the six out of eight recoupling resonances described 
by Eas. (|30l3ip that were not included in Fig. H] 



* Electronic address: 'C.Kropf@thphys.uni-heideIberg.de| 
t Electronic address: B.Fine@thphys.uni-heidelberg.de 

1 A. Redfield, Physical Review 98, 1787 (1955). 

^ A. Abragam, Principles of Nuclear Magnetism (Oxford 

University Press, Oxford, 1961). 
^ C. Slichter, Principles of magnetic resonance (Springer 

Verlag, New York, 1990). 

* J. V. Vleck, Physical Review 74 (1948). 

^ T. G. Gas, R. Griffin, and M. Levitt, J. Ghem. Phys. 89, 
692 (1988). 

® J. M. Griffiths and R. G. Griffin, Analytica Ghimica Acta 
283, 1081 (1993). 

Scholz, B. H. Meier, and M. Ernst, J. Ghem. Phys. 127, 
204504 (2007). 

* M. J. Bayro, M. Huber, R. Ramachandran, T. G. Daven- 
port, B. H. Meier, M. Ernst, and R. G. Griffin, J. Ghem. 
Phys. 130, 114506 (2009). 

^ U. Haeberlen and J. S. Waugh, Phys. Rev. 175, 453 (1968). 

° V. I. Yukalov, Phys. Rev. Lett. 75, 3000 (1995). 

^ V. I. Yukalov, Phys. Rev. B 53, 9232 (1996). 

^ See the EPAPS supplemental material for the detailed 
description of the calculations used to obtain the time- 
averaged Hamiltonians presented in this paper. 

^ D. Sakellariou, G. A. Meriles, R. W. Martin, and A. Pines, 
J. Magn. Res. Imag. 23, 295 (2005). 

* F. Bloch and A. Siegert, Phys. Rev. 57, 522 (1940). 
^ J. Winter, G. R. Acad. Sci. 241, 375 (1955). 

S. Hartmann and E. Hahn, Phys Review 128, 2042 (1962). 
^ G. P. Jaroniec, B. A. Tonnage, G. M. Rienstra, J. Herzfeld, 

and R. G. Griffin, J. Magn. Reson. 146, 132 (2000). 
^ G. P. Jaroniec, G. E. MacPhee, V. S. Bajaj, M. T. McMa- 

hon, and R. G. Griffin, Proc. Natl. Acad. Sci. 101, 711 

(2004). 



11 



=0.9 



^ = 0.1 

tin 




0.0 
-0.2 
-0.4 
-0.6 
-0.8 
-1.0 
-1.2 
-1.4 




O.Or- 

-0.2 
-0.4 
-0.6 
-0.8 
-1.0 
-1.2 
-1.4 




2 4 -4 -2 2 4 

Figure 5: (Color online) Examples of recoupling resonances for unlike spins. Plots of the resonant values of lo and the prefactors 
in the recoupled Hamiltonians pOI3ip as functions of 7j/7i for the resonant conditions indicated above the respective plots. 
All examples presented are obtained numerically for either Hi/Hq = 0.9 or Hi/Hq — 0.1. Each column of plots corresponds 
to the ratio Hi/Hq indicated above it. The plot lines consist of the actual computed points and thus in some cases have the 
appearance of dotted lines. The results for Hi/ Ho — 0.1 are qualitatively representative of the limit Hi <^ Ho- The content of 
this figure together with Fig. |4] covers all eight resonant conditions associated with Eas. (|30l3ip 



SUPPLEMENTARY MATERIAL 
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Note 1: In this supplement, we use the same notations as in the main article unless explicitly specified otherwise. 

Note 2: Some of the matrices in this supplement are computer-generated and presented as figures rather than 
equations. 



S-I. FORMALISM FOR THE CALCULATIONS OF TIME-AVERAGED SPIN-SPIN INTERACTION 

HAMILTONIANS 

We denote the full spin-spin interaction Hamiltonian in the double- rotated reference frame as 'H^(t). Our goal 
below is, first, to obtain the explicit expression for 'H^(i) and then to average it over time according to formula 

(■Hi,) = lim 7^ / H'o{t)dt. (SI) 



— r 



r— s-oo 2t 

The Hamiltonian 'H^(t) can be represented in tensor notations as follows 



JV 



i<j 

where we use the Einstein convention for summing over repeated Greek-letter matrix and vector indices, which can 
take values 1,2,3 and correspond, respectively, to the x,y,z projections in the relevant reference frame. The matrix 
^Ai!/(*i) is given by the following equation: 

A^^A^J) = P^.5ii)Cspiij)QpA:j), (S3) 

where time-independent matrix Cp,v{ij) represents the interaction coefficients of the full interaction Hamiltonian T-Ld 
in the laboratory reference frame, while time-dependent matrices Pp,v{i) and Qiiu{j) represent the transformations 
between the laboratory and the double- rotated reference frames appropriate for spins i or j, respectively. Matrix 
Cfj,i,{ij) is presented in Fig. [STl 



1-3 n^ij,^ 


3 j , X j , y 


3 ■'^i j, X ■'^i j, z 


Sym 


1-3 n^^j^y 


— 3 j, y -^i j , z 


Sym 


Sym 


1 - 3n\j,, 



Figure SI: Explicit form of the matrix Cp,y{ij). The factors riij^p = -j^r^ are time-independent material constants. The matrix 
Cfiuiij) is symmetric. 



The expressions for the transformation matrices are: 
and 
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where 



cos{ujt) sin(a;t) 0^ 
U^,^..= I -sin(t^<) cosiujt) I , (S6) 
1, 



'cos(ai) — sin(Q!i)^ 
U^..^u = I 1 

sin(Q;i) cos(Q;i) 



(S7) 



• t) - sin(ilei 

Un^..,^.^ = I sin(f7e. • t) cos(f7«,j ■ t) | . (S8) 

These are the same matrices as those defined by Eqs.(7-9) of the main article but with the spin and matrix indices 
exphcitly included. The explicit form of matrices P^ivii) and Qfj,^{j) is given in Figs. [S2]and lS3[ respectively. 





Cos [t 6J] Cos[ai] Cos[tQei] + 
Sin [t 0)] Sin [t Qgi] 


Cos[ai] Cos [t Qei] Sin[t6j] - 
Cos [t to] Sin [t Qei] 


-Cos [t Qei] Sin [tti] 


-Cos [t Qei] Sin [t OJ] + 
Cos[ta)] Cos[ai] Sin[tQei] 


Cos [t Oj] Cos [t Qei] + 
Cos[ai] Sin[ta)] Sin[tQei] 


-Sin [ai] Sin [t Qei] 


Cos [t 0)] Sin [a^] 


Sin [t w] Sin [a^] 


Cos [ai] 



Figure S2: Explicit form of the matrix P^vii) 



Cos [t w] Cos 
Sin [t 0)] Si 




h] Co 
;tQej 


s[tQej] + 
] 


-Cos [t Qej] Sin [t ijj] + 
Cos [t w] Cos[aj] Sin 


> Qej^ 




Cos [t to] Sin 






Cos [ttj] Cos [t Qej] £ 
Cos [t 10] Sin [t Qej 


in [t 0)] - 

] 


Cos [t Oj] 

Cos [aj 


Cos [t Qej] 

Sin [t 10] 


+ 

3in 


> Qej^ 




Sin [t 0)] Sin 






-Cos [t Qej] Sin 






-Sin 




Sin 


tQej 


] 


Cos [ttj] 



Figure S3: Explicit form of the matrix Q^v{j) 



The calculation of the time-average of H'jj (t) given by Eq. (jSip requires calculating 



- _lim 7^ / A^,iij)it)dt (S9) 



T~yoo 2t 



— r 



The full and time-averaged expressions for matrix A^jy{ij){t) are presented, for the case of like spins, in Section fS-III 
and, for the case of unlike spins, in Section [S-IIIl 
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To describe the different recoupling resonances we use "normalized delta-function" symbol 6(x) defined as follows: 

Six)^l' (SIO) 
otherwise. 



S-II. LIKE SPINS 



For like spins, we use the notation 

ai = aj = a (SH) 

and 

= f^ej = r!c. (S12) 



In this case, Pfj,u(i) = and, as a result, is symmetric. The time-dependent matrix A^i,{ij) is 

presented in function of the elements of C^i/(zj), Pfiuii) and Qiiv{j) in Fig. IS4I The explicit form of A^^iij) as a 
function of a, lo and Jig is given in Fig. IS5I The time-averaged value of A^i,{ij) is presented in Fig. IS6I 



Pll (CiiPn+Ci2P21 + Ci3P3l) + 

P21 (C12 Pll + C22 P21 + C23 P31) + 

P31 (Ci3 Pll + C23 P21 + C33 P31) 


P12 (C11P11 + C12P21 + C13P31) + 
P22 (C12 Pll + C22 P21 + C23 P31) + 
(Ci3 Pll + C23 P21 + C33 P31) P32 


Pl3 (CiiPii +C12P21 +C13P31) + 

P23 (C12 Pll + C22 P21 + C23 P31) + 
(Ci3 Pll + C23 P21 + C33 P31) P33 


Sym 


P12 (Cii P12 + C12 P22 + Ci3 P32) + 
P22 (C12 P12 + C22 P22 + C23 P32) + 
P32 (Ci3 P12 + C23 P22 + C33 P32) 


Pl3 (Cii P12 + C12 P22 + Ci3 P32) + 
P23 (C12 P12 + C22 P22 + C23 P32) + 
(Ci3 P12 + C23 P22 + C33 P32) P33 


Sym 


Sym 


Pl3 (Cii Pi3 + C12 P23 + Ci3 P33) + 
P23 (C12 Pl3 + C22 P23 + C23 P33) + 
P33 (Ci3 Pl3 + C23 P23 + C33 P33) 



Figure S4: Like spins : The matrix A^v{ij) in compact form. Here C^jj/ = C^i^(ij), Pin, = PiJ,v{i) and Q^jj/ = Q^i/O'). 
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Cos [a] ^ Cos [tu]^ Cos [t He] ^ 


Ci: + 


-Cos [a] Cos [t a] Cos [t !2e]^ Sin[t u] Cn + 




Cos [qJ Cos [t m]^ Cos [t Qe] Sin [a] Cn + 


2 Cos [a] Cos[t Id] 


Cos [t Qe! 


Sln[t w] Sln[t Qel Cn + 


Cos [a]^ Cos [t u]^ Cos [t Qe] Sin[t Qe] Cn - 




Cos[tw] Sin[a] Sin[tM] Sin[tQ,] Cu + 


Sln[t ai]^ Sin[t De 


]^Cn + 




Cos [t Qe] Sin[t u]^ Sin[t Q,] Cn + 




2 Cos [a] Cos [t 111] Cos[t Qe] Sin[o] Sin[t iil] C12 - 


2 Cos [a] 2 Cos [t w 


Cos [t Qe 


]^ Sln[t 111] C12 - 


Cos [a] Cos [t w] Sin[t w] Sin[t Qe]^ Cn + 




Cos[t 0)]^ Sin[a] Sin[tQ,] C12 + 


2 Cos [a] Cos[t w] 


Cos [t Qe 


] Sin[t Q,] C12 + 


Cos [a] Cos [t w] 2 Cos [t De ] ^ C12 - 




Sin[a] Sin[t OJ]^ Sin[t Q,] C12 + 


2 Cos [a] Cos[t fie 


Sln[tai; 


' Sin[t Qa] C12 - 


Cos [a] Cos [t De ] ^ Sin [t w] ^ C12 + 




Cos[a]^ Cos[t ai] Cos[tQe] C13 - 


2 Cos[t w] Sln[t M 


] Sin[t a 


= ]'Ci2- 


2 Cos [t w] Cos [t Qe] Sin[t w] Sin[t Q^] Cn + 




Cos[taj] Cos[tQe] Sin[a]^Cn + 


2 Cos [a] Cos[t w] 


Cos [t Qe! 


^ Sin [a] C13 - 


2 Cos [a] 2 Cos [t w] Cos[tQe] Sin[tw] Sin[tQ, 


Cn- 


Cos [a] Sin[t aj] Sin[t Qe] Cn + 


2 Cos [t fie] sin [a 


Sin [t bl' 


Sin[t Qe] C13 + 


Cos [a] Cos [t w] ^ Sin [t De ] ^ Cn + 




Cos [a] Cos[t Qe] Sin [a] Sin[t aj]^ C22 - 


Cos[a]^ Cos [t Qe] 


Sin[t u; 


^ C22 - 


Cos [a] Sin[t w]^ Sin[t De]^ Cn + 




Cos[ta)] Sin[a] Sinftaj] Sin[tQe] C22 + 


2 Cos [a] Cos[t w] 


Cos [t Q,; 


Sin[t w] Sin[t Q,] C22 + 


Cos [t tie]' Sin [a] Sin[t(<i] Cn - 




Cos[a]2 Cos [t Qe] Sinftai] C23 - 


Cos [t w] ^ sin [t De 


]^C22- 




2 Cos [a] Cos[tw] Cos[tQe] Sin[a] Sin[tQe] Cn - 


Cos[tQe] Sin [a] ^ Sin [t ai] C23 - 


2 Cos [a] Cos [t !2e 


2 Sln[a] 


Sin[t OJ] C23 + 


Sin[a] Sin[tid] Sin[tnel^Cn + 




Cos ' a' Cos Id] Sin [t Qe] C23 - 


2 Cos[t w] Cos [t n 


,] Sln[a' 


Sin[t C!,] C23 + 


Cos [a] Cos [t w] Cos [t fie] ^ Sin [t m] C22 - 




Cos [a] Cos[tQe] sin [a] C33 


Cos[t Q,]^ Sin [a] 


C33 




Cos[taj]^Cos[tC!,] Sin[tQe] C22 + 












Cos[a]2 Cos [t Qe] Sin[t w]^ Sin[t Qe] C22 - 












Cos[a] Cos[tid] Sin[tii)] Sin [t Q, ] ^ C22 - 












Cos [t Id] Cos [t Del ' Sin [a] C23 - 












2 Cos [a] Cos [t Del Sin [a] Sin [t bl] Sln[t 0,] C 


23 + 










Cos[tw] Sin[Q:] Sin[t Qe]^ C23 + 












Cos [t Q, ] Sin [a] ^ Sin [t tie ] C33 






Sym 






Cos [t fie]^ Sin [t u]' Cn - 




-Cos[taj] Cos[tQ(,_ Sin [a] Sin[t«] Cn + 








2Cos[a] Cos[tw] Cos[tQe] Sin[tM] SinftQ,] 


Cn + 


Cos [a] Cos[t u]^ Sin[a] Sin[t Q,] Cu + 








Cos [a] 2 Cos [t ai] ^ Sin [t Q, ] ^ Cn - 




Cos[taj]^Cos[tQe] Sin[a] Cn - 








2 Cos [t w] Cos [t Qe ] ^ sin [t w] C12 + 




Cos[tQe] Sin[a] Sin[taj]^Cn + 








2Cos[a] Cos [t w]^ Cos [t Qe] Sin[tt3e] Cn - 




2 Cos [a] Cos [t aj] Sin[a] Sin[t ai] Sin[t Q,] C12 - 








2Cos[a] Cosftfie] Sin[t ai]^ Sin[tC!o] Cn + 




Cos [a] Cos[t Qe] Sin[t aj] Cn + 








2 Cos [a] 2 Cos [t w] Sin[t OJ] Sin[t Qe]^ Cn + 




Cos[a]^Cos[ta] Sin[tQe] Cn - 








2 Cos [t Qe] Sin[a] Sin[t Ul] Sin[t Q,] Cn - 




Cos[taj] Sin[a]^ Sin[t Qe] Cn + 








2Cos[a] Cos[tw] Sin[a] Sin[tQe]^Cn + 




Cos[ta)] Cos[tQe] Sin [a] Sin[taj] C22 + 








Cos [t W]2 COS[t Qe]^ C22 + 




Cos[a] Sin[a] Sin [t a)] ^ Sin [t Qe] C22 + 








2Cos[a] Cos[t(<)] Cos[tQe] Sin[t(J] Sin[tQ,] 


C22 + 


Cos [a] Cos [t aj] Cos [t Qe] C23 + 








Cos [a] 2 Sin [t u]^ Sin [t Qe] ^ C22 - 




cos[uj bin[Eaj ^in[ti^ej C23 — 








2Cos[t!j] Cos[tQ.] Sin[a] Sin[tQe] C23 - 




Sinifl "Sin[ta] Sin[t Qe ] C23 - 








2 Cos [a] Sin [a] Sin[tw] Sin [t Q, ] ^ C23 + 




Cos [a] Sin[a] Sin[tQe] C33 








Sin [a] ^ Sin [t Q,]^ C33 






Sym 






Sym 




Cos [t ai]^ Sin [a] ' Cn + 












2 Cos [t aj] Sin[a]^ Sin[t aj] Cn + 












2 Cos [a] Cos [t aj] Sin [a] On + 












Sin[a]^ Sin[t a)]^ C22 + 












2 Cos [a] Sin [a] Sin [t ai] C23 + Cos [a] ^ C33 



Figure S5: Like spins : Explicit form of the matrix Afiv{ij) 



1 1 Cosla]' \ 

( 4 + 4 ) Cii + 

i(5(w+2ne) (2 Sin[a] - Sin[2 a] ) C13 + 

i6(W+Qe) (-l + C0S[a])2 (C11-C22) + 
/ 1 Cos [a] ^ \ ^ 1 ■ r 1 2 ^ 

(j + — C22 + J Sin[a] C33 


-i6{w+Qe) (-1 + Cos [a] ) ^ C12 + 
j(5(w+2Qe) (-2 Sin[a] + Sin[2 a] ) C23 


J- 6 (w+Qe) (-Cos [a] + Cos [a]^ - Sin [a] ^) C13 - 
i-5(2aj+Qe) (2Sin[a] -Sin[2a]) (Cn - C22) + 
i6{2a)-Qe) (2 Sin[a] + Sin[2 a] ) (Cu - C22) 


Sym 


(1^ 0031211) c^^^ 

i<5(a)+2Qe) (-2 Sin[a] + Sin[2 a] ) C13 - 

|-(5(a)+Qe) (-1 + Cos [a] )MCii - C22) + 
cosMl) C22+isin[a]2C33 


|-(5(2w+Qe) (2 Sin[a] - Sin[2 a] ) C12 + 
j6{2o)-Qg) (2Sin[a] +Sin[2a]) C12 + 
J 6 (oi+Be) (Cos [a] - Cos [a] ^ + Sin [a] ^) C23 


Sym 


Sym 


J Sin [a] ^ Cii + 1 Sin [a] ^ C22 + Cos [a] ^ C33 



Figure S6: Like spins : Time average of the matrix Afj,,^{ij). All time-averaged Hamiltonians for like spins presented in 
Section IVA of the main article follow from this expression. 
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S-III. UNLIKE SPINS 



For unlike spins, 7,; 7^ 7^. The main article associates indices i and j with two different spin species. In this section, 
in order to make formulas better readable, we change the notations as follows: 



Oil 



a 



aj = j3 



A. 



(S13) 
(S14) 
(S15) 
(S16) 



In the present case, P^u{i) 7^ Q^iv{i)i and, therefore, A^^{ij) is not symmetric. The time-dependent matrix A^^{ij) 
is presented as a function of Cpy(ij), P^_Lu{i) and Q^u{j) in Fig- IS7I The explicit form of A^^{ij) as a function of a, 
/? oj, 51 and A is given in Fig. IS81 The result of time-averaging for Af^^{ij) is presented in Fig. 



Pii (Cii Qii + C12 Q21 + Ci3 Q31) + 
P21 (C12 Qii + C22 Q21 + C23 Q31) + 
P31 (Ci3 Qii + C23 Q21 + C33 Q31) 


P12 (Cii Qii + C12 Q21 + Ci3 Q31) + 
P22 (C12 Qii + C22 Q21 + C23 Q31) + 
P32 (Ci3 Qii + C23 Q21 + C33 Q31) 


Pi3 (Cii Qii + C12 Q21 + Ci3 Q31) + 
P23 (C12 Qii + C22 Q21 + C23 Q31) + 
P33 (Ci3 Qii + C23 Q21 + C33 Q31) 


Pll (Cii Q12 + C12 Q22 + Ci3 Q32) + 
P21 (C12 Q12 + C22 Q22 + C23 Q32) + 
P31 (Ci3 Q12 + C23 Q22 + C33 Q32) 


P12 (Cii Q12 + C12 Q22 + Ci3 Q32) + 

P22 (C12 Q12 + C22 Q22 + C23 Q32) + 

P32 (Ci3 Q12 + C23 Q22 + C33 Q32) 


Pl3 (Cii Q12 + C12 Q22 + Ci3 Q32) + 

P23 (C12 Q12 + C22 Q22 + C23 Q32) + 

P33 (Ci3 Q12 + C23 Q22 + C33 Q32) 


Pll (Cii Qi3 + C12 Q23 + Ci3 Q33) + 
P21 (C12 Ql3 + C22 Q23 + C23 Q33) + 
P31 (Ci3 Ql3 + C23 Q23 + C33 Q33) 


P12 (Cii Qi3 + C12 Q23 + Ci3 Q33) + 
P22 (C12 Ql3 + C22 Q23 + C23 Q33) + 
P32 (Ci3 Ql3 + C23 Q23 + C33 Q33) 


Pl3 (Cii Qi3 + C12 Q23 + Ci3 Q33) + 
P23 (C12 Ql3 + C22 Q23 + C23 Q33) + 
P33 (Ci3 Ql3 + C23 Q23 + C33 Q33) 



Figure S7: Unlike spins : The matrix A^^{ij) in compact form. Here C^^ = Cf^viij), P/iv = PiJ,v{i) and 



Aj)- 
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Cos[al Coslffl Cos[tA] Cos[tBl^Cos[t!)l Cut 


Cos(a] Cos[(3! Cos[tal^Cos[t!)l SinftA] Cu- 


Cos|a] Cos[tfl]^Cos|tB] Sin|(3) Cut 


Cos|a] Cos[tu] Cos|ta] Si3|tA] Siiltaj Cii + 


Cos>:] Ccs|tA] Cos|ta Cos|tQ; Sin|ta] Cut 


Coslta] Sin[(3] Sin|tB] Sinita] Cut 


Cos[(5] Cos[t A] Cos[t u] Siii[t «] Sin[t !)] Cn + 


Cos|(J] Cos[tB] Sin[tA| Sin|ta| Sin|t!!l Cu- 


2Cosla] Coslta) Coslta) Sinlffl Sinjtfl] Cu - 


Siii[tAlSin[tM]^Sin[t2]Cii- 


Cos[tA] Sin[tM]'Siii[tD] Cut 


Cos[ta)^Sin[|S) SinltD) CutSin||3) Sin[ta)' SinltQ) Cut 


Cos[a] Cos[tii]^Cos[tS] Siii[tA] C12 + 


Cos[o] Cos[tA] Cos[ta]^Cos[tS] Cut 


Coslo] Cos||3) Coslta) Cos[ta) Cu- 


2Cos[al Cos[/31 Cos[tA] Cos[tM] Cos[tQ] Siii[tii] Cut 


2Cos[a] Cosfg] Cos[tM] Cos[tQ] Siii[tA] Siii[ta] Cu- 


Cosjta] Coslta] Sin|a] Sin||3) CutCosl/3) Sin[ta] Sin[tQ] Cut 


Cos[a] Cos[t()] Sia|tA] Si3|ta]^Ci2- 


Cosla] Cos[tA] Cos[t()l Sin[tal'Cu- 


Coslo] Cc5':Q] Sir|3; Sinita 'C22- 


Cos|eiCos[tA]Cos|tu]'Sin|tQlCi2- 


Cos|/i]Cos[tB]^Sin[tA]Sin[t()lCut 


Coslta] Sin 35] Sinita] Sinita] C22t 


2Cos[tu] Sin|tA] Sin|tBl Sin|tQl Cut 


2Cos[tA] Cos[tal Sin[tal Sin[t!)l Cut 


Cosia] Cos||3] Coslta; Sinita) C23- 


Cos[(5]Cos[tA]Sin[tu]^Sin[t!)]Ci2- 


Cos[|3]Sin[tA]Sin[ta]'sin[t()]Cu- 


Coslta] Sin|a] Sin|f) Sinita) C23- 


Cos[g] Cos[t A] Cos[t ii] Cos[t Q] Sinfa] C13 - 


Cos[|3] Cos[t a] Cos[t Q] Sin[a] Sin[t A] C13 - 


Cos[g] Cos|tM] Sin|tQ) C23-C0SI1S) Coslta) Sin|a] C33 


Cos[a] Cos[tA] Cos[tii] Cos[tD] Siii[/3] C13- 


Cos[a] Cos[t a] Cos[t S] Sin[S] Sin[t A] C13 1 




CosltO] Sin[a] Siii[tA] Sin[t(j] C13- 


Cos[t A] Cos[t Q] Sin[al Sin[t a] Cu - 




CoslA] Sin:S] Siilta] Si3|t2] C13 - 


Sin[ffl Sin[tA] Sin|ta Sin|ts: Cu, - 




Cosia] Cos[tu] Cos|ta] Sin[tA] Siii[tu] C22t 


Cos(a] Cos[tA] Cos|ta; Cos|ts; Sin|tB] C22t 




Cos[a] Cos[|31 Cos[tAl Cos[t01 Sin[tiii]^C22t 


Cos[a] Cos|)3j CosltQ] SinltA] Sin[tM]^C22- 




Cos[tl>i]'Siii[tA] Sin[tQ] C22- 


Cos[tA] Cos[tu]'Sin[ta] C22- 




Cos[g] Cos[tA] Cos[tlll Sin[tul Siii[tQ] C22t 


Cos[|31 Cos[tal Siii[tA] Sin(ta] Sin[tQl C22- 




Cos[tHl Cos[tB] Sin[a] Sin[tA] C23- 


Cos[t A] Cos[t u] Cos[t a] Sin[al C23 - 




Cos[S] Cos[tA] Cos[ta] Sin[al Sin[tal C23- 


Cos[ffl Cos[ta] Sin[al SinftA] Sin[tu] C23- 




Cos:nl Cos[tA] Cos[t()] Sin[g] Sin[tBl C23t 


Cosa] Cos|ta] Sin||3] SinltA] Sin[tal C23t 




Cos|tA] Cosltu] Sin[|3] Sin[tO] C23t 


Cos[tfl] Sin[(J] Sin[tA| Sin|ta] C23t 




Cos[tAlCos[tQ]Sin[olSiii[/3]C33 


Cos[tQ]Sin[alSiii[/31Sin[tA]C33 




-Cos[/3] Cos[tA] Cos[tB] Cos[t()l Sin[tu] Cu- 


-Cos[/3] Cos[tal Cos[tal Sin[tAl Sin[ta] Cut 


-Coslta] Coslta] Sinlffl Sinita] Cut 


Cos[t01 Sin[tA] Sin[tlil'Ciit 


Cos[tA] Cos[tS] Sin[tB]'Cut 


Cosia] Cos|tu]'Sinl)3] Sinita] Cut 


Cos:al Cos[/3: CosltA Cos|t u]^ Sinl S] Cut 


Cos[a] Cos|/31 Cos|tBl^Sin:tA] Sin[tc] Cu- 


Cos It a] 'Cos It a] Sin|S] Cu-Cos|ta] Sinl|3) Sinltaj'Cut 


Costal Cos[tu] Sin|tA] Sinlta; Sin|tOj Cut 


Cos a] CosltA] Cos|ta Sin|ta Sin|t5] Cut 


2Cos|a] Coslta] Sinl|3] Sinita] Sin|tQ] Cu- 


Cosiei Cos[tA] Cos|tu]^Cos|ts: Cut 


Cos;,3] Cos|ta]-Cos|ta; Sin|tA] Cu- 


CoslC] Cos|tB] Sinita] CutCos|a) Cosl/3) Cos[ta] Sin[tQ] Cu- 


2Cos[ta] Cos[tQ] Sin[tA| Sin[ta] Cu- 


2Cos|tA] Cos|ta) Cos[ta] Sin|ta] Cu- 


Coslta] Sin[a] Sinl|3) Sinita) Cut 


Cos[ffl Cos[tA] Cos[tS] Siii[t(<)l'Cu- 


CosiS] Cos|ta] SinltA] Sin|tal'Cut 


Coslta] Coslta] Sin[P] Sin|ta) C22t 


Cos;al Cos[ta]'Si3|tA: Si."|ts: Cut 


Cosli' Ccs]tA] Ccs]ta 'Sin]t5] Cu - 


Cosia] Sin||3] Sinita]^ Sinlt a] C22t 


2Ccs[a; Cos|,8] Cos;lA] Ccs[Lii] Sir.[ta] Sin[tQ] Cut 


2Ccs]a Cos],3] Coslta] Sin|tA] Sia|tB] Sin|ta] Cu - 


Cosl(3] Cosltfl] Coslta] C23tCos|a] Cos||3] Sin[ta] Sin[tB] C23- 


Cos;q] Sir.|tA] Si."|ta]-Si3|t5; Cut 


Cos ; 3] Cos ] t A] Sin ] t a] ' Sin ] 1 3] Cu t 


Sin|a] Sin)^) Sinita) Sin[ta) C23-C0SI13] Sin[a] Sin[tQ) C33 


Cos[tA] Cos|t!)] Sin[|3] Sin[tu] Cu- 


Cos[ta] Sin|(i] Sin|tA] Sin|ta] Cu- 




Cos[g] Cos[tA] Cos[tll] Siii[a] Siii[t!!] C13- 


Cos[|S] Cos|ta] Siii|a] Sin[tA] Sin[tO] C13- 




Cos[a] Cos[tA] Cos[tiil Siii[|3] Siii[tQ] C13- 


Cos[a] Cos|ta] Siii|/3] Sin[tA] Sin[tQ] Cut 




Sin[a] Sin[t A] Sin[t a] Sin[t 0] Cu - 


Cos[t A] Sin[a] Sin|t a] Sin|t a] Cu t 

Cos"tA]Cc5':a]=Cos]tQ'C::t 




Cos";i.''Ccs-:C] Sir'tA' C- t 

Cos;;3] Ccs:t.; Ccslta; Ccslta; Sinita C22- 


Cos 33' Ccsjtai' Ccsjta" Sinjt A" Sin]: a] C22t 




cos;a] ccs;t3; si3|tA; si3|ta sinits; C22t 


Coslc?: Ccs]tA: Cos]ta Sin]ta Sin]te] C22t 




Cosia] Cos[(3| CosltA) Sinita]^ Sin [tO] C22- 


Cos[a] Cos|)3| Sin|tA] Sin[tM]^Sin|tQ] C22- 




Cos[tA] Cos[tw] Cos[tO] Siii[|3] C23t 


Cos[tM] Cos[tO] Siii[|3] Sin[tA]C23- 




Cos[tal Slii[a] Siii[tA] Siii[tS] C23- 


Cos[t A] Cos[t Ii] Siii[a] Sin[t a] C23 - 




Cos[P] Cos[tA] Sin[al Sin[tal Sm\tS] C23- 


Cos[S] Sinla] Sin|tA] Sin[ta] Sin[tQ] C23- 




Cos;q] Ccs'tA' Si:.|S; Sin[:a] 3in(t01 C23t 


Cos]!' Sin]3: Sin|:A] Sinla] Sin[tQ] C23t 




Cos,tA] Sin. 3] Si.i|i3^ Sin|i5] C33 


Sin,a] Sin]/3; Sin]! A] Sin|ta] C33 




Cos[ffl Cos[tA] Cos[ta]^Sin[al Cut 


Cos[|S] Cos|ta]^Sin[(i| Sin[tA] Cu- 


Coslta)' Sin|B) Sin||3) Cut2Cos[tw] Sinfa) Sin|/3] Sin|ta] Cut 


Ccs[tii] Sin[5] SinltA Sin|ta; Cu -Cos[t a]' Sin[a] Sin[tA] Cut 


Cc.s]tA] Cos]ta] Sin]a] Sinita] CutCos]tA| Cos[ta]^ Sin[a] Cut 


Cos'6] Coslta] Sinla] Cu-Cosla] Coslta] Sin|S] Cut 


2Cos[,3] Cqs[:A] Co3":a] Sin'a] Sin[ta] Cut 


2Cos]5] Cos|r a] Sin]3] Si3]t A| Si3]ta] Cu- 


Sinia] Sinl,3; Sinlia]'C22tCosli3] Sin|a] Sinita] C23t 


Sin[a] Si3|tA] Si."|t a ' Cu -Cos|o] Cos;?] Cos[tA] Cos[tB] Cu- 


Cc5]tA] Si3]a Sin];a]^CutCos]a] Cosl(J] Cos]ta] Sin|tA] Cu- 


Cos[a) Sinl|3) Sinita] C23tCos|al Cos[|3) C33 


Cos[tA] Cos[ta] Sin[a, Sin|J] Ci3tCos[a] Sin[tA] Sin[tB] Cu- 


Cos]ta] Sia]a; Sin]i3] Sin]tA] Cu-Cos|a] CosltA] Sin|ta] Cut 




Cos[tu] Sin[a] Sin[tAl Sin[ta] C22t 


CosltA] Coslta] Sin|a] Sin[tB] C22t 




Cos[^] Cos[tA] Sm[a] Sin[ta]^C22- 


Cos[|3] Sin|a] Sin[tA] Siii|tu]'C22t 




Cos[ol Cos[tal Siii[tA] C23tCos[al Cos[g] Cos[tA] Sinfta] C23- 


Cos[a] Cos[tA] Cos[ta] C23tCos[[i| Cos[|3] Siii|tA] Siii|ta] C23- 




Cos[tA] Sin[ol Sin[S] Sin[ta] C23-Cos[a] Cos[tA] Sin[|3] C33 


Sin[a] Sinl/31 Sin[t A] Sin[tM] C23-Cos[o] Sin||3] SinltA) C33 





Figure S8: Unlike spins : Explicit form of the matrix A^„{ij) 
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-j6(ut(0+A)) ((-l + Cos[|3]) Sin[a] + (-1 1 Cos[a] ) Sin[/3] ) Cu - 


-ii5(2u-(()tA)) (-1 + Coslo]) (-l + Cos[)3]) CI2- 




ii5(wt(0-A)) ((l + Cos[/3]) Siii[o] + (-l + Cos[a]) Sin[/3]) Ci3- 


J6(2w-(0-A)) (It Cos [a]) (-ltCos[g]) C12 + 


((-1 + Cos[o]) Cos[ffl -Sin[a] Sin[3]) C13 + 


j6(w-(0-A)) ((-l + Cos[|3]) Sin[a] + (1 + Cos[a] ) Sin[|3] ) Cu - 


j6(2ut(0-A)) (-l + Cos[aj) (ltCos[g]) C12 + 


]6{2mQ,i) (-l + Cos[a]) Sin[^] (CU-C22) t 


i 5 (w- (0+A) ) ( (1 + Cos [U] ) Sin [o] + (1 + Cos [o] ) Sm [/3] ) Cu t 


i5(2u-(Q+A)) (ltCos[a]) (l + Cos[ffl) C12 + 


j6(2a)-05i) (1 + Cos[o]) Sin[3] (C11-C22) 


7i5(2u-(()tA)) (-l + Cosfal) (-l + Cos[fil) (C„-C?j) + 


j6(u+(Q+A)) ((-l + Cos[)3]) Sin[o] t (-l + Cos[a]) Sin[g]) C23- 




ii5(2u-(()-A)) (l + Cos[a]) (-ltCos[/3]) (Cii-C22) + 


j6(u+(Q-A)) ((ltCos[/3]) Sin[a] + (-1 + Cos[o] ) Sin[^] ) C23 1 




i6(24i+(Q-A)) (-l + Cos[a]) (l + Cos[S]) (Cu-C22) + 


j6(u-(Q-A)) ((-l + Cos[/3]) Siii[o] t (1 + Coslo]) Sin[/3]) 0,3- 




i6(2u-(QtA)) (l + Cos[a]) (l + Cos[ffl) (Cn-C22) 


\ 6 (u- (Q+A) ) ( (1 1 Cos [3] ) Sin [a] + (1 1 Cos [a] ) Sin [(3] ) C23 




-j6(2u-(QtA)) (-l + Cos[a]) (-l + Cos[|3]) Ci2 + 


ji5(u+(2+A)) ((-ltCos[/3]) Sin[a] + (-l + Cos[a]) Sin[3]) C13- 


-j6(2utQ,i) (-l + Cos[a]) Sin[ffl C12 1 


i6(2u-(Q-A)) (l + Cos[a]) (-ltCos[S]) Cl2- 


i6(u+(Q-A)) ((ltCos[(3]) Sin[a] + (-1 + Cos[o] ) Sin[/3]) C13 - 


i6(2u-Q,i) (l + Cos[a]) Sin[/3] C12 + 


i|5(2u+(()-A)) (-l + Cos[a]) (ltCos[/3]) Ci2 + 


\5{u-{Q-h)] ((-l + Cos[3]) Sin[o] t (1 + Cos[o] ) Sin[)3] ) C13 + 


j6(u+Qei) 


ii5(2u-(()tA)) (l + Cos[a]) (l + Cos[ffl) Ci2 + 


J 6 (w- (Q+A) ) ( (1 1 Cos [/3] ) Sin [a] + (1 1 Cos [a] ) Sin [/3] ) Cu - 


(Cos [|3] -Cos [a] Cos[/3] +Sin[ol Sin[/3]) C23 


j6(wt(0+A)) ((-l + Cos[^]) Sin[a] + (-ltCos[a]) Sin[(J]) C23 + 


j6(2u-(Q+A)) (-l + Cos[a]) (-ltCos[/3]) (C11-C22) + 




i6(wt(0-A)) ((l + Cos[|3]) Siii[o] + (-l + Cos[a]) Sin[|3]) C23- 


i6(2u-(Q-A)) (l + Cos[a]) (-ltCos[g]) (Cn-C22) + 




ji5(w-((5-A)) ((-l + Cos[|3]) Sin[a] + (1 + Cos[a] ) Sin[/3] ) C23 - 


j6(2wt(Q-A)) (-l + Cos[a]) (ltCos[/3]) (C11-C22)- 




ii5(w-((5+A)) ((l + Cos[^]) Sinfa] + (l + Cos[a]) Sin[)3]) C23 


^6(2w-((]+A)) (ltCos[a]) (l + Cos[3]) (Cn-C22) 




j6(u+QA) (Cos[a] (-l + Cosl/3]) -Sin[a] Sin[/3])Ci3 + 


-\5{2mh) (-l + Cos[)3]) Sin[a] C12 + 


\ (Sin [a] Sin[/3] Cn + 


i>5(2u+A) (-l + Cos[a])Sin[a] (Cu-C22) + 


j6(2w-A) (ltCos[/3]) Sin [a] C12 + 


Sin[a] Sin[3] C22t2Cos[aj Cos[/31 C33) 


\b(lii-ty) (l + Cos[/3]) Sin[a] (C11-C22) 


ji5((i)+QA) (Cos[a] -Cos[o] Cos[|3] +Sin[a] Sin[/i]) C23 





Figure S9: Unlike spins : time average of the matrix A^,^{ij). All time-averaged Hamiltonians for unlike spins presented in 
Section IVB of the main article follow from this expression. 



